SHIFTED SCHUR PROCESS AND ASYMPTOTICS OF 
LARGE RANDOM STRICT PLANE PARTITIONS 
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MIRJANA VULETIC* 

Abstract. In this paper we define the shifted Schur process as a 
measure on sequences of strict partitions. This process is a gen- 
eralization of the shifted Schur measure introduced in [TW] and 
[Mat] and is a shifted version of the Schur process introduced in 
[0R1| . We prove that the shifted Schur process defines a Pfafhan 
point process. We further apply this fact to compute the bulk scal- 
ing limit of the correlation functions for a measure on strict plane 
partitions which is an analog of the uniform measure on ordinary 
plane partitions. As a byproduct, we obtain a shifted analog of the 
famous MacMahon's formula. 



1. Introduction 

The basic object of this paper is the shifted Schur process that we 
introduce below. 

Consider the following two figures: 
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Both these figures represent a plane partition - an infinite matrix 
with nonnegative integer entries that form nonincreasing rows and 
columns with only finitely many nonzero entries. The second figure 
shows the plane partition as a 3-dimensional object where the height 
of the block positioned at {i,j) is given with the (i, j)th entry of our 
matrix. 

Each diagonah ((1, A;), (2, A; + 1),(3, A; + 2),...) or ((A;, 1),(A; + 
1, 2), (A; -|- 2, 3), ... ) of a plane partition is an (ordinary) partition - a 
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nonincreasing sequence of nonnegative integers with only finitely many 
nonzero elements. A strict partition is an ordinary partition with dis- 
tinct positive elements. The example given above has all diagonals 
strict and so it is a strict plane partition - a plane partition whose 
diagonals are strict partitions. 

For a plane partition vr one defines the volume Ivrl to be the sum 
of all entries of the corresponding matrix, and A{tt) to be the number 
of "white islands" - connected components of white rhombi of the 3- 
dimensional representation of tt. Then 2^^^^ is equal to the number 
of ways to color these islands with two colors. 

For a real number g, < g < 1, we define a probability measure Tig 
on strict plane partitions as follows. If tt is a strict plane partition we set 
Prob(7r) to be proportional to 2"^'^'^)g''^L The normalization constant is 
the inverse of the partition function of these weights, which is explicitly 
given by 

Proposition. (Shifted MacMahon's formula) 

oo 

TT is a strict n=l 
plane partition 

This formula has appeared very recently in [FWj and we were not 
able to locate earlier references. This is a shifted version of the famous 
MacMahon's formula: 

oo 

TT is a plane n=l 
partition 

A purely combinatorial proof of the shifted MacMahon's formula will 
appear in a subsequent publication. 

The measure described above is a special case of the shifted Schur 
process. This is a measure on sequences of strict (ordinary) partitions. 
The idea came from an analogy with the Schur process introduced in 
|0R1] that is a measure on sequences of (ordinary) partitions. The 
Schur process is a generalization of the Schur measure on partitions 
introduced earlier in [0|. The shifted Schur process we define is the 
generalization of the shifted Schur measure that was introduced and 
studied in jTW] and [MatJ. 

The Schur process and the Schur measure have been extensively stud- 
ied in recent years and they have various applications, see e.g. [B], 

[BOk], [Boi], m, m, m, m, m, m- 





^For the given example |7r| = 35 and A(7r) = 7. 
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In this paper we define the shifted Schur process and we derive the 
formulas for its correlation functions in terms of Pfaffians of a correla- 
tion kernel. We show that 971^ described above can be seen as a special 
shifted Schur process. This allows us to compute the correlation func- 
tions for this special case and further to obtain their bulk scaling limit 
when g — i> 1 (partitions become large). 

The shifted Schur process is defined using skew Schur P and Q func- 
tions. These are symmetric functions that appear in the theory of 
projective representations of the symmetric groups. 

As all other symmetric functions, Schur Px{xi,X2, ■ ■ ■) function, 
where A is a strict partition, is defined by a sequence of polynomi- 
als Px{xi,X2, . ■ ■ Xn), n eN, with a property that Px{xi,X2, ■ ■ ■ Xn, 0) = 

Px{Xi,X2, . ..Xn). 

For a partition A one defines the length 1{X) to be the number of 
nonzero elements. Let A = (Ai, A2, . . . ) be a strict partition and n > I, 
where / is the length of A. Then 



where Sn-i acts on x^+i, . . . , Xn. Schur Qx function is defined as 2 Px. 

Both {Px : A strict} and {Qx '■ A strict} are bases for Q[pi,P3,P5, . . . ], 
where pr = Yl^i power sum. A scalar product in 

Q\pi,P3,P5, • • •] is given with {Px,Q^l) = 5a,m- 

For strict partitions A = (Ai,A2,...) and /i = (/ii, /i2, . . . ), skew 
Schur functions are defined by 



where A D /i if Aj > /Xj for every i. Note that Pa/0 = -Px and Qx/is = Qx- 
This is just one of many ways to define skew Schur P and Q functions 
(see Chapter 3 of |Macj ). We give another definition in the paper that 
is more convenient for us. 

We use A to denote the algebra of symmetric functions. A specializa- 
tion of the algebra of symmetric functions is an algebra homomorphism 
A — > C. If p is a specialization and / G A then we use /(p) to denote 
the image of / under p. 

Let p = {pq,Pi,Pi, . . . ,p^) be a finite sequence of specializations. 
For two sequences of strict partitions A = (A^, A^, . . . , A"^) and fj, = 
{p^, p^, . . . , iJt^~^) we define 



W(A,p) = Qai(Po )i^AV/.i(Pi )Q\vAp't) ■ ■ - QxT /^,t-i{p1^_^)Pxt{pt). 





A D p, 
otherwise 
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Then W{X,fi) = unless 

C D /i^ C D ^ . . . D /i^-i c A^ D 0. 

The shifted Schur process is a measure that to a sequence of strict 
plane partitions A = (A^, A^, . . . , A^) assigns 

Prob(A) = ^5^W^(A,/i), 

where Z{p) is the partition function, and the sum goes over all se- 
quences of strict partitions fi = /i^, . . . , yu"^"^). 

Let X = {{xi, ti) : i = l,...,n} C N x [1, 2, . . . , T] and let A = 
(A\ A^, . . . , A"^) be a sequence of strict partitions. We say that X C A if 
Xi IS a part (nonzero element) of the partition A*' for every i = 1, . . . ,n. 
Define the correlation function of the shifted Schur process by 

p{X) = Prob(X C A). 

The first main result of this paper is that the shifted Schur process 
is a Pfaffian process, i.e. its correlation functions can be expressed as 
Pfaffians of a certain kernel. 

Theorem A. Let X C N x [1, 2, . . . T] with |X| = n. The correlation 
function has the form 

p{X) = Pf{Mx) 
where Mx is a skew-symmetric 2n x 2n matrix 

{Kxi,xjiti,tj) l<i<j<n, 
V) l<z<n<j <2n, 

where i' = 2n — i + 1 and Kx^y{ti,tj) is the coefficient of z^w'^ in the 
formal power series expansion of 

-— -J{z,ti)J{w,tj) 

2{z + w) 

in the region \z\ > \w\ if tj > tj and \z\ < \w\ if tj < tj. 
Here J{z, t) is given with 

Jit,z) = l[Fip;^;z) n ^(p™;^~'), 

t<m m<t~l 

where F{x; z) = Ylii^ + ^iz)/ (1 — Xiz). 

Our approach is similar to that of |ORlj . It relies on two tools. One 
is the Fock space associated to strict plane partitions and the other one 
is a Wick type formula that yields a Pfaffian. 



SHIFTED SCHUR PROCESS AND ASYMPTOTICS OF LARGE RANDOM SPP 5 

Theorem A can be used to obtain the correlation functions for DJlq 
that we introduced above. In that case 



Jq{t,z) 



'-q^'^z~^-q)^{q^+^nz-qU 
'-q''^z-q)oo{q''^"'z-\qU 



t > 0, 



t < 0, 



y{q'/^z-q)U-q-^+^l^z-\q), 
where 

oo 

(^;g)oo = n(i-^"^) 

n=0 

is the quantum dilogarithm function. We use the Pfaffian formula to 
study the bulk scaling limit of the correlation functions when g — 1. 
We scale the coordinates of strict plane partitions by r = logg. This 
scaling assures that the scaled volume of strict plane partitions tends 
to a constant 3 

Before giving the statement, we need to say that a strict plane par- 
tition is uniquely represented by a point configuration (subset) in 

X = G Z X Z| X > 0} 

as follows: {i, j,n{i, j)) (j — z,7r(i,j)), where is the (i,j)th 

entry of vr. 

Let 7^g, respectively 7^^, be the counterclockwise, respectively 
clockwise, arc of \z\ = R from i?e~*^ to i?e*^. 

Theorem B. Let X = {(tj, Xj) : z = 1, . . . , ra} C X be such that 
rti r, rxi — > x as r — * +0, 
ti — tj = Atij = const, Xi — Xj = Axij = const. 



a) If X > then 
where 



hm p(X) = det[ir(z,j)]-,=i, 

r— >+0 



2m J^±^ \1 + z / z^x^j+L 

where we choose if Atij > and 7^^ otherwise, where R = e"'"^'/^ 
and 

r (el"l + l)(ex-l) (el"l + l)(e>^-l) 

e = I 2el-l/2(ex + 1) ' 2el-l/2(ex + 1) " ' 

lo, otherwise. 



^The constant is equal to 7C(3)/4. 
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b) If X = and in addition to the above conditions we assume 

Xi = const 

then 

hm p(X) = Pf[M(^,j)]S=i, 
where M is a skew symmetric matrix given by 

^^11 > l<i<j<n, 



IR.B ^ ' 

1 /• (\-z\^'^i' dz 



m.3)={—\ [-r—A -J73^ \<i<n<3< 



^ iR.e ^ ' 

/■ (\-z\^'^'^' dz 



2m L± \\\Z) z-(x,,+x^,)+l 



n < i < j < 2n, 



where i' = 2n — i + l and we choose if Atij > and 7^ g otherwise, 
where R = e"''^'/^ and 6 = 7t/2. 

For the equal time configuration (points on the same vertical line) 
we get the discrete sine kernel and thus, the kernel of the theorem is an 
extension of the sine kernel. This extension has appeared in [B], but 
there it did not come from a "physical" problem. 

The theorem above allows us to obtain the limit shape of large strict 
plane partitions distributed according to Tig, but we do not prove its 
existence. The limit shape is parameterized on the domain representing 
a half of the amoe of the polynomial P{z, w) = —l+z + w + zw 
(see Section 4 for details). 

In the proof of Theorem B we use the saddle point analysis. We 
deform contours of integral that defines elements of the correlation 
kernel in a such a way that it splits into an integral that vanishes when 
q ^ 1 and another nonvanishing integral that comes as a residue. 

The paper is organized as follows. In Section 2 we introduce the 
shifted Schur process and compute its correlation function (Theorem 
A). In Section 3 we introduce the strict plane partitions and the mea- 
sure Wlq. We prove the shifted MacMahon's formula and use Theorem 
A to obtain the correlation functions for this measure. In Section 4 
we compute the scaling limit (Theorem B). Section 5 is an appendix 
where we give a summary of definitions we use and explain the Fock 
space formalism associated to strict plane partitions. 



"^The amoeba of a polynomial P{z,w) is 

{{^,io) = (log|z|,log|^z;|) e M2 I {z,w) e (C\{0})2, P{z,w) = 0}. 
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2. The shifted Schur process 

2.1. The measure. Recall that a nonincreasing sequence A = 
(Ai,A2,...) of nonnegative integers with a finite number of parts 
(nonzero elements) is called a partition. A partition is called strict 
if all parts are distinct. More information on strict partitions can be 
found in [Mac] . [Mat] . Some results from these references that we are 
going to use are summarized in Appendix to this paper. 

The shifted Schur process is a measure on a space of sequences 
of strict partitions. This measure depends on a finite sequence p = 
i.Po 1 Pi 1 ptj ■ ■ ■ ^ Pt) of specializations of the algebra A of symmetric 
functions ^. 

Let A = (A\ A^, . . . , A"^) and p = (/i^, /x^, . . . , be two sequences 

of strict partitions. Set 

W{X,fi) = Qx^{pt)P\^/^,APi)Q\vApt) ■ • •<5A^/M^-l(p^_l)A^(p^)• 
Here PA/At(p)'s and Qa//^(p)'s denote the skew Schur P and Q- 
functions, see Appendix. 

Note that W{\,fi) = unless 

C A^ D /i^ C A^ D . . . D C A^ D 0. 



Proposition 2.1. The sum of the weights W{\,n) over all sequences 
of strict partitions A = (A"*^, A^, . . . , A^) and p = {p^, /i^, . . . , p^~^) is 
equal to 

Z{p)= W H{pt,p-), (2.1) 

0<i<j<T 

where 

Hipt,pj)= E QM)Pxipj)- 

A strict 

We give two proofs of this statement. 
Proof. 1. From (15.161) and (I5.17P it follows that 

Zip) = 5^w^(A,/x) = (r+(p^)r_(p+_j---r_(p+)r+(/.r)r_(p+)^;0,^;0). 



A specialization of A is an algebra homoniorphisni A 
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We can move all F+'s to the right and all r_'s to the left using f l5.14p . 
We obtain 

z{p)= n i^(p+,P7)(r-(p?-i)---r-(Po^)r+(PT)---r+(prK^0)- 

o<i<j<r 

Then f lCTD implies (O). □ 

Proof. 2. We use Proposition 15. II The idea is the same as in Proposi- 
tion 2.1 of ^BRj . The proof goes by induction on T. 

Using the formula from Proposition 15 . 1 1 we substitute sums over A*'s 
with sums over r*~^'s. This gives 

r-i 

i=0 fJ,,T 

This is the sum of W{fi,T) with /i = {fi^ , . . . , and r = 

(r\ . . . , r^"^). Inductively, we obtain (12.11) . □ 

Definition. The shifted Schur process is a measure on the space of 
finite sequences that to A = (A^, A^, . . . , A"^) assigns 

Prob(A) = ^^iy(A,/i), 

where the sum goes over all ji = (/i^, /i^, . . . , fJp^~^)- 

2.2. Correlation functions. Our aim is to compute the correlation 
functions for the shifted Schur process. 

Let X = {{xi, U) : i = l,...,n} C N x [1, 2, . . . , T] and let A = 
(A^, A^, . . . , A"^) be a sequence of strict partitions. We will say that 
X c A if part of the partition A*' for every i = 1, . . . ,n. 

Definition. Let X C N x [1,2, .. .T]. Then the correlation function 
of the shifted Schur process corresponding to X is 

p{X) = Prob(X C A). 

We are going to show that the shifted Schur process is a Pfaffian pro- 
cess in the sense of |BR] . that is, its correlation functions are Pfaffians 
of submatrices of a fixed matrix called the correlation kernel. This is 
stated in the following theorem. 

Theorem 2.2. Let X C N x [1,2, .. . T] with \X\ = n. The correlation 
function is given with 

p(X) = Pf (Mx) 
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where Mx is a skew- symmetric 2n x 2n matrix 

{K^^^^^{ti,tj) l<i<j<n, 
{-lp'K,^^_,,^{ti,ty) l<t<n<j <2n, 

{-IY^'+-.'K.,^,^.^^, (ti,, V) n<i<j<2n, 

where i' = 2n — i + 1, j' = 2n — j + 1 and Kx,y{ti, tj) is the coefficient 
of z^w^ in the formal power series expansion of 

K{{z,ti),{w,tj)) := J2 K^^y{ti,tj)z''wy = 

in the region \z\ > \w\ if ti > tj and \z\ < \w\ if ti < tj. 
Here J{z, t) is given with 

Jit,z) = l[Fip-;z) n F{pi;^-'), (2.2) 

t<m m<t—l 

where F is defined with ( f5. 3\} . 

Proof. The proof consists of two parts. In the first part we express the 
correlation function via the operators r+ and r_ (see Appendix). In 
the second part we use a Wick type formula to obtain the Pfaffian. 
Let = Po = and let to = 0. 

First, we assume I < ti < t2 < ■ ■ ■ < tn < T . Using formulas fl5.16p . 
f l5.17p and (15.111) we get that the correlation function is 



^(-A\Yl'^-(p^)'^+(p^)Yl{^'^^^^*x^ n r-(p+)r+(p„))t;0,t;0 



where the products of the operators should be read from right to left 
in the increasing time order. 

Thus, the correlation function is equal to nr=i('~-'-)^' times the co- 
efficient of nr=i in the formal power series 



Z{p) 



n r-(p;t.)r+(p-)n(2^("^)^(^^) n ^~(p^)^+(p 

m 

m=ti-i 



where if) is given by (I5.12p . 

We use formulas (I5.14p and ( 15.15^ to put all r_'s on the left and 
r+'s on the right side of HILi '^4'{ui)ilj{vi). Since T^vq = vq, see (I5.13p . 
we obtain 

Y[{J{ti,Ui)J{ti,Vi))/Y[2ij{ui)ij{vi)v9,V9y (2.3) 

1=1 \ i=l ' 
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Thus, ii 1 < ti < t2 < ■ ■ ■ < tn ^ T then the correlation function 
p(X) is equal to nr=i(~-'^)^' times the coefficient of YYi=i'^i*'^i~^' 
the formal power series (12.31) . 

Now, more generally, let n E Sn such that 1 < < t^(2) < ■ ■ ■ < 
tn{n) < T, then the correlation function p{X) is equal to nr=i(~-'-)^' 
times the coefficient of nr=i in the formal power series 

Y[{J{ti,Ui)J{ti,Vi))( Jj2V'(M^(i))V'(t^,r(i))^^0,^^0 V (2.4) 
i=l ^ i=l ' 

The above inner product is computed in Lemma 12.41 and Lemma 12.51 
that will be proved below. Then by Lemma 12.51 we obtain that in the 
region |'U7r(„)| > \vn{n) \ > ■ ■ ■ > 1^77(1) I > expression (12.41) is equal 

to Pi{A) with 

{K{{ui, U), {uj, tj)) I <i < ] <n, 
K{{u,,U),{v'^,t'.)) l<i<n<j<2n, 
K{{v[,Q,{v'j,t'^)) n<t<j<2n, 

where i' = 2n — i + 1, j' = 2n — j + 1 and K is as above. 

Let (2/1, . . . ,?/2n) = (a;i, . . . . . . , -a;i) and (zi, . . . , 2;2„) 

,fi). By the definition of the Pfaffian, 



pf(Mx) =n(-i)'''E^g^(«)n^^«. 



^Vjs 1 

i=l 



where the sum is taken over all permutations a 

' 1 2 ■ ■ ■ 2n - 1 2n . , ^ • ■ . ■ ■ r 

such that ti < ■ ■ ■ < tn and ts < js for 



«i Ji ••• ^n 3 
every s. 

Also, 



Pf(A) = ^sgn(a)n^(^..,%) 



2n 



Finally, since 



we get 



yi,---,y2n&^ a s i=l 



p{X) = \[{-m\[uTvr-.Vi{A)] 

i=l 1=1 

p(X) =Pf(Mx). 



□ 
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Lemma 2.3. 

2n 
i=2 

(2.5) 

Equivalently, 

where A is a skew symmetric matrix with Aij = {tpkj'ipkiV^, v%) for i < j ■ 

Proof. Throughout the proof we use (■■■^0,^0) = (■ ■ ■ ) to shorten the 
notation. 

We can simphfy the proof if we use the fact that the operators ipk 
and ipl add, respectively remove e^, but the proof we are going to give 
will apply to a more general case, namely, we consider any i/j^^s that 
satisfy 

iJkV9 = 0, fc<0, (2.6) 
^poV9 = aovii, (2.7) 
rk = h^-k, (2.8) 

i^ki'i + i'lipk = CkSk-i, (2.9) 

for some constants oq, &fc and Ck- These properties are satisfied for ipk 
we are using (see Appendix). 

The properties (I2.6p - fl2.9l) immediately imply that 

iJki^k = 0, k^Q (2.10) 



and 

because for k ^ —I 



{^kiJi)=0, k^-l (2.11) 



ii^ki^i) = { -bi{ipkV0,ip-iv,i) I > 0, 
ao{ilJkV<!),v,i) 1 = 0. 

Observe that if /ci, . . . , k2n-i are all different from then 

{^k,„.,---^k,) = 0. (2.12) 

This is true because there is > such that {ki, . . . , k2n-i) has 
odd number of elements whose absolute value is k. Let those be 
fcj^, . . . , kj2^_-^. If there exist kj- and kj.^-^ both equal either to k or —k 
then ijJkir,-! ■ ■ - ^ki = by (12. 9p and (12.101) . Otherwise, {kj^, . . . , kj.2^_J 
is either {k, —k, . . . , —k, k) or (—A;, k, . . . ,k, —k). In the first case we 
move ipk to the left and use (12.61) and (12.81) and in the second we move 
ip-k to the right and use (12. 6p . 



12 



M. VULETIC 



Also, observe that if {kn, . . . ,ki) contains t zeros and if those appear 
on places jt, ■ ■ ■ ,ji then 

{^k^ ■ --^k,) = ■ ■ -4,, ■ ■ -4,, ■ ■ ■ (2.13) 

This follows from ( 12. 7p and ( 12. 9p by moving ipk^^ to the mth place from 
the right. 

Now, we proceed to the proof of the lemma. We use induction on n. 
We need to show that L{k2n, ■ ■ ■ ,ki) = R{k2n, ■ ■ ■ ,ki) where L and R 
are the left, respectively right hand side of (12.51) . Obviously it is true 
for n = 1. We show it is true for n. 

If A;i < then L = R = hy ([21]). 

If fci = then let t be the number of zeros in k2n, ■ ■ ■ ,ki and let 
those appear on places jt, . . . , ji = 1. Then 

t 



m=l 



while 

t 

R = ^i-'^y'{i^k2r,---'i'k,^---'^k2){lpk,^M 



\=2 
t 



1=2 m=2 

t t 

m=l i=2 

If t is even then L = R. If t is odd then 2n — t is odd then 
(^fen ■ ■ -4,, ■ ■ ■4,, ■ --M = by dm and thus L = R=0. 

Finally, we can assume ki > 0. If ki ^ —ki for every i G {2, . . . ,n} 
then 

L{k2n, ...,ki)= R{k2n, ...,ki)=0. 

This is true because i? = by (12.111) and L = since it is possible to 
move ipki to the left and then use (12.61) and (12. 8p . 

So, we need to show that L = R if ki > and there exists i G 
{2, ...,«} such that ki = —ki. 

First, we assume ^2 = — fci. Then 

L = (V^fe2„---V'fe3)(V'fc2V'fel), 

because ip^kipk^di = {ip-k'ipk)v(D. On the other hand 
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because for every i > 2 we have that {'>pk2n ' ' ' i^ki ■ ■ ■ ipk2) = by (12. 6p . 
Thus, L = R. 

To conclude the proof we show that 

L{k2ny • • • 5 ks, kg-i, . . . , ki) = R{k2n^ ■ • • 1 ksy kg-i, . . . , fci) 

^ (2.14) 
L{k2ny • • • ) ks_i, kg, . . . , ki) = R{k2n-! ■ ■ ■ ) ks_i, kg, . . . , ki). 

Then we use this together with the fact that there is i such that ki = 
—ki and that L = i? if i = 2 to prove the induction step. 

L{k2n, ks-i, ks,..., ki) = (^fc2„ ■ ■ ■ ■ ■ V'fci) = 

= -L{k2n,---,ks,ks-i,...,ki) + Ckjks-ks-i{i^k2n ■ ■ ■ '4^ks'^ks-i ■■■V'fci) 
R{k2ny • • • ; kg-i, kg, ■ ■ ■ , ki) = 

2n 
i^SjS— 1 

= —R{k2n-, • • • ; ^s; • • • 5 ^l) 

2n 

+Cfc,4.,-fc._i X] (-l)*(^fc2n--4r-4.4.-i--V^fci)(^fc,V'fci)- 

i=2 

i^s, s—1 

Since, by the inductive hypothesis 

2n 

j=2 
tT^SjS— 1 

we conclude that (I2.14p holds. 



□ 



Lemma 2.4. Let n G Sn- 

JJV'(M7r(i))V'K(i))^^0,f0 ) = Pf(*7r), 

where '^j^ is a skew- symmetric 2n-matrix given with 

{'ilj^{{ui,i),{uj,j)) l<i<j<n, 
{v'pj')) l<t<n<j <2n, 
MWi,^'), n<i<j <2n, 
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where i' = 2n — i + 1, j' = 2n — j + 1 and 



Proof. First we show that the statement is true for tt = id. From 
Lemma [231 we have 



i=l 



UiVi 



2rx n 

i=2 i=l 
n 

i=l 

n 



1=2 



Then by the expansion formula for the Pfafiian we get that 

{i){Un)lp{Vn) ■ ■ ■ 1p{Ui)ll){Vi)v%, V%) = Pf (v4(u„, W„, ...,Ui, Vi)), 

where A{un, Vn, ■ ■ ■ ,Ui,Vi) is a skew symmetric 2n x 2n matrix 



a{Un,Vn) 

-a{Un, Vn) 

-a{un,ui) -a{vn,ui) 

_ -a{un,vi) -a{vn,vi) 



a{vn,ui) a{vn,vi) 

a(ui,t>i) 
—a{ui,vi) 



with 



a{z,w) = {i>{z)i>{w)vq,,V(i^). 

Thus, the columns and rows of A appear "in order" 
Mji, . . . , Ml, fi. Rearrange these columns and columns "in or- 
der" Ml, • • • , M„, . . . , V2, Ml. Let B be the new matrix. Since the 
number of switches we have to make to do this rearrangement is equal 
to n{n — 1) and that is even, we have that Pf(A) = Pf(-B) and 

-a(M2,Mi) ... -a(M2,Mi) a(Mi,Mi) 

a(M2,Mi) ... a(M2,M2) a(M2, Mi) 

B 

a(M2,Mi) -a(M2,M2) ... a(M2, Mi) 

_ -a(Mi,Mi) -a(M2,Mi) ... -a(M2,Mi) 

This shows that B = "^id. 
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Now, we want to show that the statement holds for any it E Sn- We 
have just shown that 

{ij{Unytp{Vn) ■ ■ ■ 'tp{ui)ij{vi)V(l,, V^) = Pf . . . , Ui, Vi)), 

where B is given above, i.e. 

{B{{ui,i),{ujJ)) l<i<j<n, 
B{{u,,t),iv'^,f)) l<t<n<j <2n, 
B{{vi,i'),{v'j,j')) n<i<j<2n, 

where i' — 2n — i + 1, f — 2n — j + 1 and 



B{{z,i),{w,j)) 



{ij{z)ij{w)v(i,,VQ) i>j, 

-{ljj{w)ll!{z)V(l,,V(A) i<j. 



Then 

^7r(n)) ^7r(n)) ■ ■ ■ ) ^7r(l)) 

Change the order of rows and columns in S(it^(„), i'7r(n), ■ ■ ■ , it7r(i), 'i'7r(i)) 
in such a way that the rows and columns of the new matrix appear 
"in order" -Ui, ^2, • • • , Vn, . . . ,V2,Vi. Let C be that new matrix. The 
number of switches we are making is even because we can first change 
the order to v„, . . . , ui, Vi by permuting pairs {ujVj), and the num- 
ber of switches from this order to Ui, U2, . . . ,Un, . . . , ^2, t'l is n{n — l) 
as noted above. Thus Pf(C) = Pi{B). Then C = because 

(B{{ui,7r-\i)),{uj,7r-\j))) 1 < i < j < n, 
C{i,j) = I B{{ui,7r-\i)), {v'j,7r-\j'))) l<i<n<j<2n, 
[B((vl,7r-\i')),(v'j,7r~\j'))) n< i < j < 2n. 

□ 

Lemma 2.5. Let tt e Sn. In the domain |M7r(n)| > \'V-K{n)\ > ••• > 

|''^7r(l)| > |'^7r(l)| 



{Yl'^{u^{{))'^(v7r{j))v$,Vg,) = Pf(*), 



i=l 



where "if is a skew- symmetric 2n-matrix given with 

{'ijj{ui,Uj) 1 <i < j <n, 
ip{ui,v'j) l<i<n<j<2n, 
i^iv'i,v'j) n<i<j<2n, 
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where i' = 2n — i + 1, j' = 2n — j + 1 and 

z — w 



A{z + w)' 



Proof. This is a direct corollary of Lemma 12.41 and a formula given in 
Appendix: 

z — w 

{'4){z)i){w)vq,,V(ii) = — — ■ for |z| > 

A{z + w) 

It is enough to check for 1 < i < j < n, because other cases are 
similar. If vr^^(i) > 7r^^(j), respectively vr^^(i) < 7r^^(j) then 

(mj, j)) = i){ui,Uj). (2.15) 

in the region \ui\ > respectively \ui\ < \uj\. This means that (12.151) 
holds for the given region |M7r(n)| > \'Vn{n) \ > ■ ■ ■ > \Utt(i)\ > |'^^7r(i)|- D 

Remark. Equation (12.51) appears in |Matj . Lemmas 12.41 and 12.51 are not 
of interest in |Matj for n ^ id because in the case of the shifted Schur 
measure one does not need to consider the time order. 



3. Measure on strict plane partitions 

In this section we introduce strict plane partitions and a measure on 
them. This measure can be obtained as a special case of the shifted 
Schur process by a suitable choice of specializations of the algebra of 
symmetric functions. Then the correlation functions for this measure 
can be obtained as a corollary of Theorem 12.21 Using this result we 
compute the (bulk) asymptotics of the correlation functions as the par- 
titions become large. One of the results we get along the way is an 
analog of MacMahon's formula for the strict plane partitions. 

3.1. Strict plane partitions. A plane partition tt can be viewed in 
different ways. One way is to fix a Young diagram, the support of 
the plane partition, and then to associate a positive integer to each 
box in the diagram such that integers form nonincreasing rows and 
columns. Thus, a plane partition is a diagram with row and col- 
umn nonincreasing integers associated to its boxes. It can also be 
viewed as a finite two-sided sequence of ordinary partitions, since each 
diagonal in the support diagram represents a partition. We write 
TT = (A~"^^, . . . , A"^, A°, A^, . . . , A^^), where the partition A° corresponds 
to the main diagonal and A'^ corresponds to the diagonal that is shifted 
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I I I I I 
I I I I I I 



X X X X X X X X 



Figure 1. A strict plane partition vr = 

((2), (3, 2), (4, 3, 1), (5, 3, 2), (3, 2), (2, 1), (1), (1)) 

by ±k, see Figure [U Every such two-sided sequence of partitions rep- 
resents a plane partition if and only if 

(3.1) 

skew diagrams A'^^"'^ — A*^ (A*^ — A*''^^) are horizontal strips. 

A strict plane partition is a finite two-sided sequence of strict parti- 
tions with property (13. ip . In other words, it is a plane partition whose 
diagonals are strict partitions. As an example, a strict plane partition 
corresponding to vr = ((2), (3, 2), (4, 3, 1), (5, 3, 2), (3, 2), (2, 1), (1), (1)) 
is given in Figure [TJ 

We give one more representation of a plane partition as a 3- 
dimensional diagram that is defined as a collection of 1 x 1 x 1 boxes 
packed in a 3-dimensional corner in such a way that the column heights 
are given by the filling numbers of the plane partition. A 3-dimensional 
diagram corresponding to the example above is shown in Figure |2l 

The number |7r| is the norm of vr and is equal to the sum of the filling 
numbers. If vr is seen as a 3-dimensional diagram then Ivrl is its volume. 

We introduce a number A{tt) that we call the alternation of vr as 

tr+i -n 

i=l i=0 

where A"^^"*"^ = A^-'"^^^ = 0, and / and a are as defined in Appendix, 
namely /(A) is the number of (nonzero) parts of A and a(A — fi) is the 
number of connected components of the shifted skew diagram A — /i. 

This number is equal to the number of white islands (formed by white 
rhombi) of the 3-dimensional diagram of the strict plane partition. For 
the given example A{tt) is 7 (see Figure [3]). In other words, this number 
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< 



r 





is equal to the number of connected components of a plane partition, 
where a connected component consists of boxes associated to a same 
integer that are connected by sides of the boxes (i.e. it is a border 
strip filled with a same integer). For the given example there are two 
connected components associated to the number 2 (see Figure [3]). 



of the connected components. We state this fact as a proposition and 
prove it. 

Proposition 3.1. Let n be a strict plane partition. Then the alternate 
A{7t) defined with liS. 2\) is equal to the number of connected components 
ofn. 

Proof. We show this inductively. Denote the last nonzero part in the 
last row of the support of vr by x. Denote a new plane partition obtained 
by removing the box containing x with n'. 






defined by 03.21) is equal to the number 
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We want to show that A{tt) and A{7c') satisfy the relation that does 
not depend whether we choose fl3.2p or the number of connected com- 
ponents for the definition of the alternate. 

We divide the problem in four cases I, II, III and IV shown in Figure 
m and we further divide these cases in several new ones. The cases 
depend on the position and the value of xl and xr. 





Figure 4. Cases I, II, III and IV 

Then using (13.21) we get 

A{7r') = y4(7r)+contribution from x^^+contribution from x/j+change of 

Let us explain this formula for case I when xl = xr = x in more 
detail. 

Let Xxj^, Xx and A^^^ be the diagonal partitions of vr containing xl, 
X and Xr, respectively. Let be a partition obtained from A^, by 
removing x. Then 

contribution from xl = a(A^ — A^^^^) — a{Xx — Xx^) = 0, 
contribution from xr = a{Xxj^ — A^.) — a(A^^ — A^^) = 1, 
change of /(A°) = /(A°(7r)) - l{X%n')) = 0. 
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For all the cases the numbers are 



A{tt') = A{tt) + 



I 


II 


III 


IV 




0+1+0 


0+0+1 


0+0+0 


0+0+0 


XL = X,Xr = X 


-1+0+0 


-1-1+1 


1+0+0 


0-1+0 


Xl > X,Xji > X 


0+0+0 


0-1+1 






Xl = X,Xji > X 


-1+1+0 


-1+0+1 






XL > X,Xr = X 



It is easy to verify that we get the same value for A{7i') in terms of 
AItt) using the connected component definition. Then inductively this 
gives a proof that the two definitions are the same. □ 

We can give a generalization of Lemma ISTTl Instead of starting from 
a diagram of a partition we can start from a type of a diagram shown 
in Figure [5] (connected skew Young diagram) and fill the boxes of this 
diagram with row and column nonincreasing integers such that integers 
on the same diagonal are distinct. Call this object a skew strict plane 




Figure 5. SSPP 

partition (SSPP). We define connected components of SSPP in the 
same way as for the strict plane partitions. Then we can give an analog 
of Lemma [3?T] for SSPP. We will not use this result further in the paper. 

Lemma 3.2. Let ir be a SSPP 

TT = (■ • ■ C A*-2 D ■ ■ ■ D A*-i C ■ ■ ■ C A*° D ■ ■ ■ D A*i C • • ■ C D ■ ■ ■ )• 

Then the alternate A{7i) defined by 

oo oo oo 

A{7c)= J2 «(A\A-i)- ^(^*")+ E ^(^*""')' 

i=— oo i=— oo i=— oo 

with 

I a{fi — A) ^ D X 
is equal to the number of connected components o/vr. 
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Proof. The proof is similar to that of Lemma I3.1[ □ 

3.2. The measure. For q such that < g < 1 we define a measure on 
strict plane partitions by 

Pi^) = (3.3) 

where Z is the partition function, i.e Z = This mea- 

sure has a natural interpretation since 2^'^'^^ is equal to the number of 
colorings of the connected components of vr. 

This measure can be obtained as a special shifted Schur process for 
an appropriate choice of the specializations p. Then Proposition 12.11 
can be used to obtain the partition function: 



Proposition 3.3. 



IT is a strict n=l 
plane partition 



This is the generating formula for the strict plane partitions and we 
call it the shifted MacMahon's formula since it can be viewed as an 
analog of MacMahon's generating formula for the plane partitions that 
says 



vr is a plane n=l 
partition 



Before we show that this measure is a special shifted Schur process let 
us recall some facts that can be found in Chapter 3 of |Macj . We need 
the values of skew Schur P and Q functions for some specializations 
of the algebra of symmetric functions. They can also be computed 
directly using (15.11) . 

If p is a specialization of A where Xi = s, X2 = = . . . = then 



Qx/^,{p) 



2a(x-^l) ^\x\-\^l\ A D p, A - /i is a horizontal strip, 

otherwise, 

2a(A-M)-z{A)+z{^)^|A|-||.| A D /i, A - p is a horizontal strip, 

otherwise, 



where, as before, a(A — /i) is the number of connected components of 
the shifted skew diagram A — /i. In particular, if p is a specialization 
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given with xi = X2 = ■ ■ ■ = then 



1 ^ = fi, 

otherwise, 

1 A = 

otherwise. 



In order to obtain the measure (13.31) as a special shifted Schur process 
we set 

p+ : Xi = g-(2n+i)/2^ = = . . . = n < —1, 
p- : xi = X2 = . . . = n<-l, , . 

p-:xi = g(2"+i)/2, = X3 = . . . = n > 0, ^ ^ ^ 

: Xi = X2 = . . . = n > 0. 



This measure is supported on strict plane partitions viewed 
as a two-sided sequence (. . . , A~", . . . , A°, . . . , A", . . . ). Indeed, for 
any two sequences A = (. . . , A~", . . . , A°, . . . , A", . . . ) and p = 
(. . . , . . . , . . . , p", . . . ) we have that the weight is given with 

00 

n=— 00 

where only finitely many terms contribute. Then W{\,p) = unless 

• ■ ■ c A-" C ■ ■ • C A° D ■ ■ ■ D A" D ■ • • , 

skew diagrams A"+^ — A" for n < and 
A" — A"+^ for n > are horizontal strips, 

and in that case 



;,a(A"-A"-i) (-2n+l)(lA"|-lA"-i|)/2 



A" n<0, 

Xn+l n > 0, 



n=— 00 
00 



2a(A"~^-A")-/(A"-i)+/(A")^(2n-l)(|A"-i|-|A"|)/2 



'q 

n=l 



Thus, the given choice of p's define a shifted Schur process (or a limit of 
shifted Schur processes as we explain in a remark below) that is indeed 
equal to the measure on the strict plane partitions given with (13. 3p . 
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Proposition 12 . II allows us to obtain the shifted MacMahon's formula. 
If p+ is Xi = s, X2 = = . . . = and p~ is Xi = t, X2 = X3 = . . . = 
then 



H{p\p-) = l[- 



1 + XiVj 



_l + st 

x=p+,y=p~ ^ 



Thus, for the given specializations of p^'s and p~'s we have 

1 + g 1 + g2 1 + g3 



Z{p) = \[H{pt.p- 



I — q \ — (f 1 — 

1 + 1 + 

1 — g2 1 — 

1 + g^ 1 + g^ 



n 

n=l 



g3 1 - g4 

1 + g"^ " 
1 - g" 



Remark. A shifted Schur process depends on finitely many specializa- 
tions. For that reason, measure (13. 3p is a limit of measures defined as 
shifted Schur processes rather than a shifted Schur process itself. For 
every T let specializations p^ be as in (13.41) if |n| < T and zero other- 
wise. They define a shifted Schur process whose support is St that is 
the set of strict plane partitions with A" = for every \n\ > T. The 
partition function for this measure is Yl-KeSr and is bounded 

by IKf^f^^)"- Let S be the set of all strict plane partitions. Then 
the Tth partial sum (sum of all terms that involve g*" for m < T) of 
^^gg 2'^('^)gl'^l , which is equal to the Tth partial sum of X^TreSr S'^'^'^-'g''^' , 
is bounded by rKi^f^^)"- Hence, Yln^s'^'^^^^^^'^^ converges. Therefore, 
E^eSr 2^^^^^'^' ~^ E^es^^^^^g'''' as T ^ 00. Thus, the correlation 
function of the measure (13. 3p is the limit of the correlation functions 
of the approximating shifted Schur processes as T — > cxd. 

Our next goal is to find the correlation function for the measure ( 13. 3p . 
For that we need to restate Theorem 12.21 for the given specializations. 
In particular, we need to determine J(t, z). When p is such that xi = 
s,X2 = X3 = ■ ■ ■ = then (see (15. 3p ) 

F{p] z) = . 

1 — sz 
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Thus, for our specializations (see (12.21) ) 

1 + 



J{t,z) 



m>t 



m>0 



1 + ^-+1/2^-1 (_^i/2^-i. q)^(qt+i/2z. q)^ 



t > 



q 



771+1/2^-1 



nf 

m>0 



1 + 

771+1/2^ 



m>—t 



i-q'/'z;qUiq-'^y'z-';q), 

1 + g^'+V^^-l (gl/2^; g)^(_g-m/2^-l. 
g7n,+l/22; — 1 



t < 0, 



(3.5) 



where 



{z;q)^ = l[il-q-z) 



77=0 



is the quantum dilogarithm function. 

It is convenient to represent a strict plane partition vr = 
(. . . , A~^, A°, A^, . . .) as a subset of 

X = {(t,x) eZxZ\x>0}, 

where {t, x) belongs to this subset if and only if a; is a part of A*. We 
call this subset the plane diagram of the strict plane partition vr. 

Corollary 3.4. For a set X = {(tj, Xi) : i = 1, . . . ,n} C X representing 
a plane diagram, the correlation function is given with 

p{X) = Pf(Mr) (3.6) 

where is a skew-symmetric 2n x 2n matrix and Mx^{i,j) is given 
with 



■ I (ti ; tj 



l<i<n<j< 2n, 



(3.7) 



(-l)"''+"^'^-x„-.y (t7',t,0 n<z<j<2n, 
where i' = 2n — i + 1 and K{{ti, x), {tj, y)) is the coefficient of z^w^ in 



the formal power series expansion of 

z — w 



K{{t,,z),{t„w)) 



2(z + w^ 



■J(ti, z)J(tj,w) 



in the region \z\ > \w\ if t^ > tj and \z\ < \w\ if ti < tj. Here J{t,z) is 
given with hS. 5\) . 
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4. ASYMPTOTICS OF LARGE RANDOM STRICT PLANE PARTITIONS 

In this section we compute the bulk hmit for shifted plane diagrams 
with a distribution proportional to 2^^'^'^q^^\ where vr is the correspond- 
ing strict plane partition. In order to determine the correct scaling we 
first consider the asymptotic behavior of the volume |7r| of our random 
strict plane partitions. 

4.1. Asymptotics of the volume. The scaling we are going to 
choose when computing the limit of the correlation functions will be 
r = — log q for all directions. One reason for that lies in the fact that 
r^|7r| converges in probability to a constant. Thus, the scaling assures 
that the volume tends to a constant. Our argument is similar to that 
of Lemma 2 of [QRl] . 

Proposition 4.1. If the probability of n is given with /13. 3\) then 
The convergence is in probability. 

Proof. We recall that if E{Xn) c and Var(X„) then X„ — c in 
probability. Thus, it is enough to show that 

and 

Var{r^\n\) ^0, r ^ +0. 
First, we observe that 

and 

Var{\.\) = E{\nf)-Ei\n\f = ' / = ^^d-D- 

Since Z = nn>i(l5^)"' ^^'^^ 



m 



- Zj 



^ n>l ^ ^ ' m>l I ± I m>l ^ 
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Then 

m>l ^ m>l,fe>0 A;>0 ^ ^ 

because 



m>l ^ 



Now, 



^ (1 - g2fe+i)3 - /^(2/,-+l) + o(/-(2/,- + l))V ^ (2A; + 1)3 ^ 



r 

Then by the uniform convergence 

hm r=^£;(|7r|) = 2 V hm r''L^±J-^L^ = 2 V — 



Finally, since 



it follows that 



ifc>0 ^ 



^(^'kl) ^ ^C(3), r^+O. 
For the variance we have 

since by L'Hopital's rule 

!c(3) ^ ita £p = to M. 

Thus, 

V ar {r^ \7r\) ^0, r ^ +0. 

□ 
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4.2. Bulk scaling limit of the correlation functions. We compute 
the limit of the correlation function for a set of points that when scaled 
by r = — logg tend to a fixed point in the plane as r +0. Namely, 
we compute the limit of (13.71) for 

rti — > r, rxi — > x as r — > +0, 

ti — tj = Atij = const, Xi — Xj = Axij = const, 
where x ^ 0. 

In Theorem 14.21 we show that in the limit the Pfaffian (13. 6p turns 
into a determinant whenever x > and it remains a Pfaffian on the 
boundary x = 0. 

Throughout this paper 7^_g ilufi) stands for the counterclockwise 
(clockwise) oriented arc on \z\ = R from i?e~*^ to i?e*^. 

More generally, if 7 is a curve parameterized by i?(0)e"^ for (j) G 
[— 7r,7r] then 7^" (7^) stands for the counterclockwise (clockwise) ori- 
ented arc on 7 from R{6)e~^^ to R{6)e^^ . 

Recall that our phase space is X = {{t,x) G Z x Z | a; > 0} . 

Theorem 4.2. Let X = {{ti, Xi) : i = 1, . . . ,n} G X be such that 
rti T, rXi — > X as r ^ +0, 
ti — tj = Atij = const, Xi — Xj = Axij = const, 
a) If X > then 

jim^p(X) = det[K{i,j] 

where 



KUj) = —[ (^—^] "—l—dz, 

where we choose 7^^ if Atij > and 'j'^g otherwise, where R = e"'"^'/^ 
and 

f (e'"l + l)(ex-l) (el"l + l)(e^-l) 

larccos 2eM/2(ex + 1) ' 2eM/2(ex + 1) " ' (4.1) 

\^0, otherwise, 

b) If X = ^ '^^^ iiT' addition to the above conditions we assume 

Xi = const 

then 

limp(X) = Pf[M(z,j)]S=i, 
1 — >+o 
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where M is a skew symmetric matrix given by 



2ni 
1 



z 



dz 



l + Z^ 
Z ^ 



2711 



l + z 

1 



dz 



l + Z 



At,-/,-/ 



dz 



I < i < j < n, 
l<i<n<j< 
n < i < j < 2n, 



where i' 



2n — i + l and we choose if ^tij > and'y]^g otherwise, 
where R = e"''^'/^ and 6 = 7r/2. 

This implies that for an equal time configuration (points on the same 
vertical line) we get 

Corollary 4.3. For X = {{t, Xi),i = 1, . . . ,n} C X such that 
rt r, rXi ^ x as r +0, 

3j ^ X j — i^j — CO TX S 

where x > 0, 

/,^N -, rsin(6'Axjo 
hm p(X) = det — 4 ^ 

r-^+O I TcAXij 

where 6 is given with ( [^.jp . 

Remark. The kernel of Theorem 14.21 can be viewed as an extension of 
the discrete sine kernel of Corollary 14.31 This is one of the extensions 
constructed in Section 4 of [B] , but it is the first time that this extension 
appears in a "physical" problem. 

The limit of the 1-point correlation function gives a density for the 
points of the plane diagram of strict plane partitions. We state this as 
a corollary. 

Corollary 4.4. The limiting density of the point (r, x) of the plane 
diagram of strict plane partitions is given with 

Q 

piT,x)= , lim K3Diit,x),{t,x)) 



lim 

rt—fT,rx^Xt 



71 



where 9 is given with ( [y^.ip . 

Remark. Using Corollary 14.41 we can determine the hypothetical limit 
shape of a typical 3-dimensional diagram. This comes from the obser- 
vation that for a strict plane partition 



x{t,x) 



pij, s)ds, 
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if r > and 



y{'r,x) 



p{t, s)ds + r, 



X, 



p(r, s)ds — T, 



p{t, s)ds, 



z{r,x) = X 

if r < 0. Indeed, for a point (x, y, z) of a 3-dimensional diagram, 
where z = z{x,y) is the height of the column based at {x,y), the 
corresponding point in the plane diagram is (r, x) = {y ~ 2;, z). Also, 
for the right (respectively left) part of a 3-dimensional diagram the 
coordinate x (respectively y) is given by the number of points of the 
plane diagram above (r, x) (i-e. all (r, s) with s > x) ^-nd that is in 
the limit equal to p{t, s)ds. 



The hypothetical limit shape is shown in Figure 6 (a 





(a) The limit shape 



(b) The amoeba of -1+z+w+zw 



IS 



Figure 6. 

Recall that the amoeba of a polynomial P(z,w) where {z,w) G 

{(e,cu) = (log|^|,log|«;|) gM^ I {z,w) G (C\{0})^ P{z,w)=0}. 

The limit shape of the shifted Schur process can be parameterized 
with {^,Lj) = (r/2,x/2) for {^,uj) G S) where 2D is the shaded region 
given in Figure 6(b) The boundaries of this region are u = 0, u = 
log[(e« + l)/(e« - 1)] for ^ > and u = log[{e~^ + l)/(e-« - 1)] for 
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^ < 0. This region is the half of the amoeba oi —1 + z + w + zw for 
u = log \w\ > 0. 

Proof. (Theorem 14.21) Because of the symmetry it is enough to consider 
only r > 0. 

In order to compute the limit of fl3.7p we need to consider three cases 

1) = -^Xi,X2('^l) ^2) 

when rxi —>■ x, rti r, r ^ +0, q = and Ati^2 and Axi^2 are fixed. 

Since we are interested in the limits of 1), 2) and 3) we can assume 
that rti > e^/^, for every i. 

We start with 2). By the definition, 

= // -— -J{ti,z)J{t2,w)——^-—dzdw, 

|z| = (l±e)e^/2 
|io|={l^e)e'^/2 

where we take the upper signs if ti > t2 and the lower signs otherwise. 
Here, we choose e G (0, 1 — g^/^) since J(t, z) is equal to its power 
series expansion in the region g^/^ < \z\ < g^*~^/^. With the change of 
variables w t— > —w we get 



where 

z + w J{ti,z) / z \~-^^ / W 



(27ri)2 J J 2zw{z-w)J{t2,w) \e^/^J 

|^| = (l±e)e^/2 
|«,|=(l^e)e^/2 



dzdw. 



We consider cases 1) and 3) together because terms in the Pfaffian 
contain and K in pairs. 

Using the definition and a simple change of coordinates w ^ —w we 
have that 

where 

1 ff z + w J(ti,z) f z / w \-^2 

dzdw 



{2Tiif JJ 2zw{z -w) J{t2,w) \e^/^J \e^/^ 
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and 

1 ff z + w J{h,z) / z y-i ( ^ ^ ^ 
(27rz)2 J J 2zw{z-w)J{h,w)\e-/V \e-/V ' 

|^| = (l±e)e^/2 
|^| = (l=f:e)e^/2 

where the choice of signs is as before: in the first integral we choose 
the upper signs if ti > t2 and the lower signs otherwise; similarly in 
the second integral with ^3 and t^. 

In the first step we will prove the following claims. 



Claim 1. 



1 r 1 /I - e-^;z\^*''' / z \-Axi,2 



where we choose the plus sign if ti > t2 and minus otherwise and 6 is 
given with (14. ip . 

Claim 2. 

Ito = lim -L(_l)..«.« f 1^24 (4.) rf,. 

r-.+o ^^+o27rz^ ' zJ{ti,z)\e^/^) 

where we choose the minus sign if ti > t2 and plus otherwise, where 
6 is given with (14. ip . (We later show that the limit in the right hand 
side always exists.) 



Claim 3. 



lim/__= hm^ / '-^^(A.r"^'dz, 



r^+o 27rz L± zj(t4^,z) \e^/'^ 

where we choose the plus sign if ^3 > ^4 and minus otherwise, where 
9 is given with (14. ip . (We later show that the limit in the right hand 
side always exists.) 

In the next step we will show that Claims 2 and 3 imply that the 

limits of J++ and / vanish when r +0 unless X = 0. We state this 

in two claims. 

Claim 4. If % > then limj-^+o I++ = and if x = then 

1 [ 1 /I - e^^z\^*''' / z \-(^i+^2) 



72,^/2 



i^o^++ 2m J z\l + e-^z) Ve^/^. 



where we pick 7^/2 ^^1,2 > and otherwise. 
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Claim 5. If x > then lim,r^+o I = and if x = then 



where we pick 7^,, ,^ if t^t-xA > and 7 ,^ otherwise. 

We postpone the proof of these claims and proceed with the proof 
of Theorem 14.21 

a) If X > then the part of the Pfaffian coming from ++ and 

blocks (two n X n blocks on the main diagonal) will not contribute to 
the limit. This is because every term in the Pfaffian contains equally 

many and K factors. Now, by (14.21) we have that ■ K'^t = 

const (—1)'^2+^^ ■ /++/ and then by Claim 4 and Claim 5 we have that 

This means that the Pfaffian reduces to the determinant of H — block, 
A 

' det A. By Claim 1 we have that 



because Pf 



-A 



2mJ^± \l + z, 



Ati,2 



and it is easily verified that e'^*-'"-' prefactors cancel out in the determi- 
nant. 

Thus, if X > then for 

rti T, —>■ X as r ^ +0, 
ti — tj = Atij = const, Xi — Xj = Axij = const, 
we have that 

p(X)^det[ir(^,j)]-,=i, 
where K{i,j) is given in the statement of the theorem. 

b) Now, if X = and Xi does not depend on r, then by Claims 1, 4 
and 5 we have that 

= (^lY2^-r{x^+X2)J_ f (^—^\ i dz 

^ ^ 27ri L± \l + z J 2^1+^2+1 ' 

^ ' 27ri J^± z\l + e-^z 
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= ('_X)^3e'^^^^"^''*^— — / I - - I dz 

^ ^ 2m J ^± \1 + Zj ;2-fe+^4) + l 

and 

;^l,2 _ ^~r{x,-X2)J_ [ n-z \'^'''' 1 

+ 27ii J^± \l + zj ^(^i-:j;2)+i 

As before, it is easily verified that e^*^"' ^ prefactors cancel out in the 
Pfaffian. 

Hence, if x = then for 

rti —>■ r, rXi — > x as r — >■ +0, 

tj — tj = Atij = const, Xi — Xj = Axij = const, Xj = const 
we have that 

p(X)^Pf[M(z,j)]S=i, 

where M{i,j) is given in the statement of the theorem. 
It remains to prove Claims 1 through 5. 

Proof. (Claim 1) In order to compute the limit of /+_ we focus on its 
exponentially large term. Since [z + w)/2zw{z — w) remains bounded 
away from z = w, z = and w = 0, the exponentially large term comes 
from 

J{ti,z) / z / w exp [log J(ti, — xi(log2 — r/2)] 



J{t2,w) \e^/'^ / exp [log J(t2, w^) — a^2(log w — t/2)] 

(4.3) 

To determine a behavior of this term we need to know the asymptotics 
of log(2;; g)oo when r — > +0. Recall that the dilogarithm function is 
defined by 



(1 — zY 
1 ^ 

n=l 

with the analytic continuation given by 



dilog(l - .) = / ^rft 



with the negative real branch cut. Then 

log(z;g)oo = --dilog(l - z) + 0(dist(l, r^lO < r < 1)"^), r +0. 
r 

(see e.g. [B]). Hence, (14.31) when r — > +0 behaves like 
exp ^ [S{z, T, x) - S{w, r, x)] + 0(1), 
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S{z, r, x) = -dilog(l + e^'^z) - dilog(l ) + dilog(l - e^'^z) 



+dilog(l + 



Xilogz- r/2). 



The real part of S{z, r, x) vanishes for 1^1 = e"^/^. We want to find the 
way it changes when we move from the circle \z\ = e"^/^. For that we 
need to find ^ 

—ReS{z,T, x)- 

On the circle, ReS'(z, r, x) = implies that the derivative in the tangent 
direction vanishes, i.e. 



y- 



dReS dReS 



X- 



0. 



dx dy 

The Cauchy-Riemann equations on \z\ = e^l"^ then yield 



^^Re5(z,r,x) = 
Simple calculation gives 

45(.,r,x) 

which implies that 
d 



c/ReS* (iReS* 
X — ; h y 



dx 



-X + log 



^is{z,r,x). 



dy 

[l + e-^z){z + l) 
[l-e-^z){z-l) 



z—S{z,T,x) 
dz 



'X + log 



z + 1 



Then 



d 



z^S{z,T,x) > iff < 
dz 



z-1 

z+1 



for \z\ 



for \z\ 



r/2 



(4.4) 



radius -^zj where R = e^l'^ . If 



z-1 



One easily computes that e^l'^ 

/2. li 

+ l)(e>^- 1) 



is a circle with center and 



> 1 



(4.5) 



2e^/2(ex + 1) 

then this circle does not intersect \z\ = e^l'^ and ^Re5'(2;, r, x) is neg- 
ative for every point z such that \z\ = e^^^. Otherwise, this circle 
intersects \z\ = e^^^ at z = e'^^^e^*^ where 9 is given with fl4.ll) . Thus, 
the sign of ^ReS'(2;, r, x) changes at these points. 

Let 72 and 7^; be as in Figure [71 We pick them in such a way that 
the real parts of S{z, r, x) and S{w, r, x) are negative everywhere on 
these contours except at the two points with the argument equal to 
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±9 (if (14.51) is satisfied tliey are negative at tliese points too). This is 
possible by (IQ]) . 




Figure 7. Contours 7^ and 7^, 

Tlie reason we are introducing 7^ and 7^0 is that we want to deform 
the contours in /+_ to these two and then use that 

f S{z) 

/ e~dz -^0, r +0. 

J-y 

whenever S has a negative real part for all but finitely many points on 
7 . When deforming contours we will need to include the contribution 

coming from the poles and that way the integral will be a sum 

of integrals where the first one vanishes as r — +0 and the other 
nonvanishing one comes as a residue. 

Let f{z, w) be the integrand in /_|__. We consider two cases a) ti > ^2 
and b) ti < ^2 separately. 

Case a) ti > t2- We omit the integrand f{z,w)dzdw in the formulas 
below. 
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(27ri)2 
1 

(27ri)2 
h- 

Ii + 



/ / — 27ri / Res^=;j/(2;, tu)^^; 



2z J{tuz 



2; / 2; 



=r/2 



^2; 



zJ{t2,z) Ve 



-Aa;i,2 



dz. 



The first integral denoted with 7i vanishes as r — > +0, while the 
other one goes to the integral from our claim since 



lim 



r-*+0j{t2,z) \l + e-'^Z 

Case b) ti < t2 is handled similarly 



|z|=(l-e)e^/2 |«,|=(l+e)e^/2 Iz |«,|=(l+e)e^/2 



(27ri)^ 




-2m 



Res^=2/(2;, w)dz 



^ I2 + 



2ni 



2m 



2z J{ti,z) 
-2z^ J{t2,z) Ve 

Zj{t2,z) 



□ 



Proof. (Claim 2) The exponentially large term of /++ comes from 



J{ti,z) 

J{t2,w) 



^1 / w \-;^2 exp [log J(ti, — Xi(log2; — r/2)] 



/ z / w \ 



exp [log J(t2, w) + X2(logw - t/2)] 



that when r — > +0 behaves hke 



where 
S(z,T, x) 



exp-[S{z,r,x) +T{w,T,x)] + 0{l), 
r 



-dilog(l + e^'z) - dilog(l ) + dilog(l - e^'z) 

z 
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+dilog(l + -) -x(log^-r/2) 

and 

T(w, r, X) = dilogfl + e^'^w) + dilogfl ) — dilogfl — e~'^w) 

w 

-dilog(l + -) - xi^ogw - r/2). 
w 

Real parts of S{z, r, x) and T{w, r, x) vanish for 1^1 = e"^/^ and \w\ = 
e"^/^, respectively. To see the way they change when we move from 
these circles we need -^ReS{z,T,x) and -j-ReT{w,T,x)- For \z\ = 
e^/^ and |w| = e'^^'^ they are equal to ■^-^S{z,t,x) and ^^T(w,r, x), 
respectively. 

The needed estimate for S is (14.41) above. 

Similarly, 

w-^T{w, r, x) > iff < 
dw 

Hence, if (14. 5 p holds then both -^ReS{z, r, x) and ^ReT(z, r, x) are 
negative for every point z such that \z\ = e"^/^. Otherwise, the sign of 
■j-ReS{z, T, x) changes a.t z = e^l'^e^'^^ ^ while the sign of ^ReT(w, r, x) 
changes at w = e'^/^e^*^'^"^-'. 

We deform the contours in /++ to 7^ and 7^ shown in Figure [8] 
because the real parts of S and T are negative on these contours (except 
for finitely many points). 







w 



1 



w + 1 



for \w\ = e^^"^ . 




Figure 8. Contours 7^ and 7, 
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As before, we distinguish two cases a) ti > t2 and b) ti < t2- 
For a) ^1 > t2 (where we again omit f{z,w)dzdw) 



1 



/ / 



1 



(27ri)2 J J (27ri)2 

|z|=(l+e)e^/2 |«)|=(l-e)e-^/2 Iz |u)|=(l-e)e-^/2 



(27ri)2 

_ 1 
(27ri)2 

_ 1 
(27ri)2 

^'~2^/ 



+ 



27ri / ReSt(,=^;/(2;, w)^^; 




7^ •'7io 



— 27ri / Res^=^/(2;, 



e-r/2,,r-e 

2^; J{ti,z) f z / z \-^2 



-2^2 J(^2,^) 



For the same reason as in the proof of Claim 1 we have that /i — > 
when r — > +0. Thus, hm^-v+o I++ — hm^-^+o h, where 

TTiL ^ J(t2,;^) Ve-/^; Ae-Z^J 



27ri 

_ 1 

= (-1) 

For b) ti < t2 
1 



1 J(ti,-2;) 



^ -zJ{t2,-z) \e^/^ J Ve-/2 
6^/2, e 



—2; 



-2; 



-12 



dz 



X1+X2+ 



2" A-,,, 



1 J{t2, Z) ( Z \-{xi+X2) 



zJ{ti,z) 



dz. 



/ / 



(27ri)2 J J (27ri)2 

|z|=(l-e)e^/2 |w| = (l+e)e^/2 7^ j«)| = (l+e)e^/2 



1 



(27ri)2 




2m I ReSu,=zf{z,w)dz 
22; J(ti,^) / 2; / 2; \-^2 



e' / ,7r- 



-2^2 
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Thus, linir-^+o -^++ = where 

1 J{ti, z) f z 



h = 



21:1 J - Zj(t2,z) 

27ri —z J{t2, —z) 



=r/2 



or/2 



-X2 



dz 



(-1) 



X1+X2 + 1 . j 



1 J(t2, z) ( Z N -(2:1+3:2) 



zJ{ti,z) 



Therefore, we conclude that 

hm /. 



++ 



hm /, 



with 



where we choose j~r/2 5, if > ^2 and 7^/2 g otherwise. □ 

Proof. (Claim 3) The exponentially large term of / comes from 

J(t3, z) (^Y^ /_u^Y^ ^ exp [log J(t3, z) + a;3(log z - r/2)] 
J{U,w)\e'l'^) exp[log J(i4,w) -X4(logw -t/2)] 

that when r — > +0 behaves hke 



where 
S{z,T, x) 

and 



exp - [S{z, r, x) + T{w, r, x)\ + 0(1), 



-dilog(l + e"^z) - dilog(l ) + dilog(l - e"" z) 

z 



+dilog(l + -) + x(log2;-r/2) 



T{w,T,x) = dilog(l + e'^^iy) + dilog(l ) - dilog(l - e~''t(;) 

w 

-dilog(l + -) + xi'^ogw - r/2). 
w 

As before, real parts of S{z., r, x) and T{w., r, x) vanish for \z\ — e"^/^ 
and |w| = e"^/^, respectively. 



Since 



z^>5(z,T,x) > iff e^/2> 



2 - 1 



z+ 1 



for \z\ — e^^"^ 
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and 

w^T{w, r, x) > iff > 
aw 

we deform contours in / to 7^ and 7^ shown in Figure M 



w + 1 



w 



for \w\ = e"^/^, 



Y 

t w 



Figure 9. Contours 7^ and 7^ 



Using the same reasoning as in the proof of Claim 1 and 2 we get 
that 



hm 



hm J, 



with 



1 J(t3, -2) / ;z 



27ri L± zJ(tA,z) Ve"^/^ 



-0:4 



where we choose 7^/2 ^ if ^3 > ^4 and « otherwise. 



□ 



Proof. (Claim 4) We start from the integral on the right hand side in 
Claim 2. 



Because 



we focus on 



hm — = 

r^+O J{z,t2) 

Z \ -{^1+^2) 



1 - e-^z 



Ati 



1 + e-^z 
exp 



(4.6) 



where 



S{z,T,x) = -2xlogz + xr. 
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Assume x > 0- Then ReS'(z, r, x) > if and only if \z\ < e"^/^. 

For Ati 2 > we deform the contour 7~^/2 g to 7^ (see Figure [T0|) . 
Using the same argument as in the proof of Claim 1 we get that 
lim^^+o = for Ati^2 > 0. A similar argument can be given 
for Ati,2 < 0. 




Figure 10. Contour 7^ 

For X = the claim follows directly from Claim 2 changing z —z 

and using (14. 6p . □ 

Proof. (Claim 5) The proof is similar to the proof of Claim 4. □ 

The proof of Theorem 14. 21 is now complete. □ 



5. Appendix 

5.1. Strict partitions. A strict partition is a sequence of strictly de- 
creasing integers such that only finitely many of them are nonzero. 
The nonzero integers are called parts. Let throughout this section 
A = (Ai, A2, . . . ) and /x = (/xi, /i2, • • • ) be strict partitions. 

The empty partition is = (0, 0, . . . ). 

The weight of A is |A| = Yl ^i- 

The length of A is /(A) = # of parts of A. 

The diagram of A is the set of points (i, j) G such that 1 < j < A^. 
The points of the diagram are represented by 1 x 1 boxes. The shifted 
diagram of A is the set of points E 71? such that i < j < Aj + i — 1. 
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The diagram and the shifted diagram of A = (5, 3, 2) are shown in 
Figure [Hi 



Figure 11. Diagram and shifted diagram of (5,3,2) 



A partition is a subset of A if /ij < Aj for every i. We write jj C X 
or A D /i in that case. This means that the diagram of A contains the 
diagram of /i. 

If A D /i, the (shifted) skew diagram A — is the set theoretic differ- 
ence of the (shifted) diagrams of A and fi. 

If A = (5, 3, 2) and fi = (4, 1) then A D /i with the skew and the 
shifted skew diagram A — /i as in Figure [121 



r 



I r 



Figure 12. Skew diagram and shifted skew diagram of 
(5, 3, 2) -(4,1) 

A skew diagram is a horizontal strip if it does not contain more than 
one box in each column. The given example is not a horizontal strip. 

A connected part of a (shifted) skew diagram that contains no 2 x 2 
block of boxes is called a border strip. The skew diagram of the example 
above has two border strips. 

If A D /i and A — is a horizontal strip then we define a(A — fi) as the 
number of integers i > 1 such that the skew diagram X — fi has a box 
in the ith column but not in the (z + l)st column or, equivalently, as 
the number of mutually disconnected border strips of the shifted skew 
diagram X — fi. 

Let P be a totally ordered set 

P = {1 < 1' < 2 < 2' < • ■■}. 

We distinguish elements in P as unmarked and marked, the latter being 
one with a prime. We use |j9| for the unmarked number corresponding 
to p e P. 
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A marked shifted (skew) tableau is a shifted (skew) diagram filled 
with row and column nonincreasing elements from P such that any 
given unmarked element occurs at most once in each column whereas 
any marked element occurs at most once in each row. Examples of a 
marked shifted tableau and a marked shifted skew tableau are 
5 3' 2' 1' 1 1' 1 

3 2' V 2' V 

11 11 

For X = {xi, X2, ■ ■ ■ , Xn) and a marked shifted (skew) tableau T we 
use X to denote where is equal to the number of 

elements p in T such that \p\ = i. 

The skew Schur function Q\/^ is a symmetric function defined as 

Q\/t,{Xl,X2,...,Xn) = {^'^^ ' (5.1) 

lU, otherwise, 

where the sum is taken over all marked skew shifted tableaux of shape 
X — fi filled with elements p E P such that \p\ < n. Skew Schur 
Pa//, function is defined as Pa/^ = 2'(^)-'(^)Qa//.- Also, one denotes 
Pa = Pa/0 and Qx = Qa/o- 

This is just one of many ways of defining Schur P and Q functions 
(see Chapter 3 of |Macj ). 

We set 

1 + XiUj 



H{x,y) = ll^^ (5.2) 
and 

oo 

1 + XiZ 



i=l * 

F{x;z) is denoted with Qx{z) in [Mac]. 

Let A be the algebra of symmetric functions. A specialization of 
A is an algebra homomorphism A ^ C. If p is a specialization of A 
then we write Pa/^(p) and Qx/^{p) for the images of Pa/^ and Qx/^, 
respectively. Every map : {xi,X2, ■ ■ ■) — > (fli, 02, . . . ) where Oj G C 
and only finitely many a^'s are nonzero defines a specialization. For 
that case the definition (15.11) is convenient for determining Qx/^i{px)- 
We use H{px,Py) and F^p^, z) for the images of (15.21) and (15. 3p under 
Px ® Py and px, respectively. 

We recall some facts that can be found in Chapter 3 of [Mac] : 



A strict i,j 



H{x, y)= Yl QA(x)PA(y) = n (5-4) 
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Qx{x,z)= Qx/^{x)Q^{z), (5.5) 

strict 

Px{x,z)= PxU^)P,i^), (5.6) 

fi strict 
u strict 

P,ix)P,ix) = Y /4^a(x), (5.8) 



A strict 

fA 



for /^^, G Z. 
Proposition 5.1. 

A strict T strict 

Proof. Let if = H{x,y)H{x,u)H{z,y)H{z,u). Then by (15.41) we have 

= 5^QA(x,^)PA(y,t^). 
A 

We can compute H in two different ways using (15. 5p . (15. 6p . (15.71) and 
(15.81) . One way we get 

H = 5^gA(x,^)PA(l/,w)= $^gA/M(a;)QM(^)^A/.(l/)P.(w) 
A A,/i,v 



E 



X1^Vm(2^)^a/z.(z/) 

A 

The other way we get 

H = H{x,y)YQAx)PMQpi^)Ppiy)Qr{^)Pr{u 



Q^{z)P,{u). (5.9) 



}i{x,y)YQA^)Pp{v){Qp{mr{z)){PMPM) 

H{x, y) Y QAx)P,{y) Y 2-'^'^)2'('')2'W/^^,g,(;.) Y nrPM 

<T,p,T fl V 

H{x,y) Y QM(^)P.(^)($^2'(^)2'('')2-'(^)/;,p,(y))(5^/;,g.(x)) 

H{x,y) Y Qp{z)PAu)Pp/r{y)Qu/r{x) 
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E 



H{x,y)^P^/r{y)Qu/ri 



X] 



Q^{z)P,{u). (5.10) 



Now, the proposition follows from (15.91) and (I5.10p . □ 
We give another proof of this proposition at the end of Section 15.21 

5.2. Fock space associated with strict partitions. We introduce 
a Fock space associated with strict partitions. We follow [Mat] . See 
also ID.TKMj . 

Let V be a vector space generated by vectors 

vx = A ex^ A ■ ■ ■ A exi, 

where A = (Ai, A2, . . . , A;) is a strict partition. In particular, = 1 
and is called the vacuum. 

For every G Nq = NU{0} we define two (creation and annihilation) 
operators ipk and ip^. These two operators are adjoint to each other for 
the inner product defined by 

{vx,v,) = 2-'^^^5x,,- 

For k = 0, let 

(_1)HA) 
-ipovx = ^Ivx = Vx. 

For k > 1, the operator ipk adds Cfc, on the left while ipl removes e^, 
on the left and divides by 2. More precisely, 

'4'kVx = Cfc A Vx, 

^fc^A = 2 '^'^'^'^^1 A • ■ ■ A A ■ ■ ■ Ca, . 

1=1 

These operators satisfy the following anti-commutation relations 

^|J^^pJ+Mi = 0, (^,j) eN2\{(0,0)}, 

iij:^* + rj^: = o, (^,j)gn2\{(o,o)}. 

Observe that for any i G N 

jvx/2 if^GA, 
i^i-^iVx = < ^ . (5.11) 

I otherwise. 

Let ip{z) be a generating function for ipk and ipk- 

^{z) = Y,^kz', (5.12) 
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where 
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Then the following is true 
and 



(-l)V-fc k<-l. 



1)72 A; = -/ > 1, 
k = l = 0, 
otherwise, 




k=0 
Z — W 



w 



for \z\ > \w\ 



A{z + w)' 

For every odd positive integer n we introduce an operator and its 
adjoint q;_„ = a*: 

an = y^(-l)^^fc-n^-fc, 



kez 



a. 



-n = ^(-l)Vn-fe^fc- 



fcez 



Then the following commutation relations are true 



n 



for all odd integers n and m, where [a, f3] = aj3 — j3a. 

We introduce two operators and its adjoint F_. They are defined 

by 

2pn(a;) 



E 

.n=l,3,5,. 



n 



F-|-(a;) = exp 
Here p„ G A are the power sums. Then one has: 

F+(x)f0 = 



T+{y)T_{x) = H{x,y)T_{x)T+{y), 



(5.13) 
(5.14) 



V±{x)i){z) = F{x; z^^)ij{z)T±{x), (5.15) 

where H{x,y) and F{x; z) are given with (15. 2p and (15. 3p . respectively. 

The connection between the described Fock space and skew Schur P 
and Q functions comes from the following: 
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T_{x)Vf, = Y Qx/f,ix)vx, (5.16) 

A strict 

T+{x)vx = ^ Px/f,{x)Vf,. (5.17) 



fi strict 

We conclude this appendix with another proof of Proposition 15. 1[ 
Proof. (Proposition 15. ip 

{T_{x)Vf„T_{y)vx) = ( '5A/M(a;)^^A, Q^/^^v)^^ 

A strict A strict 



A strict 



^ T strict T strict 



T strict 



By (15.141) we have 

{T+{y)T_{x)v^,v^) = H{x,y){T_{x)T+{y)Vf,v^) 
which imphes 

Y Q^/t^i^)P^/Ay) = H{x,y) Y Qy/r{x)P^/r{y) 

A strict T strict 

□ 
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